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Abstract 



Generic slow-fast systems with only one (time-scaling) parameter 
on the two-torus have attracting canard cycles for arbitrary small val- 
ues of this parameter. This is in drastic contrast with the planar case, 
where canards usually occur in two-parametric families. In present 
work, general case of nonconvex slow curve with several fold points is 
considered. The number of canard cycles in such systems can be effec- 
tively computed and is no more than the number of fold points. This 
estimate is sharp for every system from some explicitly constructed 
open set. 
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For the planar case (i.e. (x, y) G M 2 ), there is a rather simple description of its 
behavior for small e. It consists of interchanging phases of slow motion along 
stable parts of the slow curve M := {(x,y) \ f(x,y,0) = 0} and fast jumps 
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1 Introduction 



Consider a generic slow-fast system: 
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along straight lines y = const. (See e.g. [5].) Given additional parameters, 
depending on e, one can observe more complicated behavior: appearance of 
duck (or canard) solutions (particularly limit cycles), i.e. solutions, whose 
phase curves contain an arc of length bounded away from uniformly in e, 
that keeps close to the unstable part of the slow curve [4]. . 

In [1], Yu. S. Ilyashenko and J. Guckenheimer discovered a new kind of 
behavior of slow-fast systems on the two-torus. It was shown that for some 
particluar family with no auxiliary parameters there exists a sequence of 
intervals accumulating at 0, such that for any e from these intervals, the 
system has exactly two limit cycles, both of which are canards, where one 
is stable and the other unstable. Yu. S. Ilyashenko and J. Guckenheimer 
conjectured that there exists an open domain in the space of slow-fast systems 
on the two-torus with similar properties. Here we prove this conjecture, and 
provide almost complete description for bifurcations of canard cycles on the 
two-torus. In particular, we give sharp estimate for the number of canard 
cycles in such systems. 

Our main results are the following ones. Consider the system (1) and 
assume that the phase space is the two-torus: 

(x,y) e T 2 = M 2 /(27rZ 2 ). (2) 

Assume that the speed of the slow motion is bounded away from zero (g > 0), 
the slow curve M is a smooth connected curve, and its lift to the cover- 
ing coordinate plane is contained in the interior of the fundamental square 
{\x\ < it, \y\ < 7r}. We also assume that all fold points of the slow curve 
(i.e. the points of M where the tangent line to M is parallel to x-axis) are 
nondegenerate (i.e. the tangency rate is quadratic). In this case, the number 
of fold points is finite and even: let us denote it by 2N. 

Theorem 1.1. For any generic slow-fast system on the two-torus with the 
properties described above, under some additional nondegenericity assump- 
tons, the following properties hold. There exists a positive number k < N and 
a sequence of intervals accumulating to zero, such that for every e belonging 
to one of these intervals the system has exactly k attracting and k repelling 
limit cycles. All these cycles are canards, and make exactly one turn along 
the y-axis during the period. The measure of their basins of attraction or re- 
pulsion is bounded away from uniformly in e. Finally, for any sufficiently 
small e > 0, the number of limit cycles that make one turn along the y-axis 
does not exceed 2k. 

Theorem 1.2. There exists an open set in the space of slow-fast systems 
on the two-torus for which the maximal number k of pairs of canard cycles 
reaches its maximal possible value N. 
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The paper is organized as follows. In section 2 we provide an heuristic 
description of the phenomena discussed. In section 3 some preliminary results 
about the Poincare map are stated. Section 4 gives an overview of the proof of 
Theorem 1.1. This proof relies on some auxiliary results, that are discussed 
in sections 5 and 6. Section 7 is devoted to construction of system with 
maximal number of canards and proof of Theorem 1.2. 

Due to size limitations, we omit technical details from presented proofs. 
We refer the reader to the works [2, 3] for more detailed discussion. 

2 Description of the phenomena 

In this section we provide heuristic description of the phenomena discovered 
by Ilyashenko and Guckenheimer. In what follows, we will assume that x- 
axis of fast motion is vertical, and y-axis is horizontal. The slow motion is 
directed from the left to the right. 

We consider first the simplest case: M is a convex curve and therefore it 
has exactly two fold points (i.e. N = l). 1 The right one is called jump point 
and the left one is reverse jump point. Consider a strip B in the phase space 
that contains M and bounded by vertical circles that pass through the fold 
point (see Fig. 1). We call it the base strip. In more generic (nonconvex) 
case the base strip is defined as the minimal vertical strip that contains M. 

Fix some vertical cross-section T = {y = const} that does not interset M. 
We will assume without loss of generality that T = {y = — 7rJ = {y = n}. 
Consider some point w ^ M from the interior of the base strip B. Trajectory, 
passing through this point, in forward time is attracted quickly to the stable 
part of the slow curve, then moves slowly to the right until reaches the jump 
point, then "jumps" and continues slow motion along the y-axis, making 
about 1/e rotations along the x-axis before it intersects T (call this phase 
after-jump rotations). For given e, denote the point of the first intersection 
with T by R(e). 

In backward time, the trajectory is quickly attracted to the repelling part 
of the slow curve, then moves slowly to the left until reaches the reverse jump 
point, jumps, and continues slow motion along the y-axis while rotating along 
the x-axis, up to the intersection with T. Denote the point of the intersection 
by L(e). This trajectory is canard, because it has a segment which is close 
to the repelling part of the slow curve. 

As e > decreases, "fast" parts of the trajectory become more verti- 
cal, and the number of rotations during the after-jump motion increases. 

^dn fact, only the latter condition matters: any system with N = 1 can be considered 
as a system with convex slow curve. 
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Figure 1: Canard solution of the system with convex slow curve: the base 
strip is shadowed 



Therefore, the point R(e) moves upwards, and L(e) moves downwards. By 
continuity, there exists E\ such that for e = E\ these two points coincide: 
R{e\) = L(e\). This gives us canard limit cycle. As e > continues de- 
creasing, new coincidense occurs for some e — 62, < e% < Si, and so on. 
Therefore, for the sequence of parameter values e = accumulating to 0, 
the system has canard limit cycles. By choosing initial point close enough to 
the reverse jump point, it is possible to make these cycles stable. When we 
perturb initial point slightly, corresponding values of Ek also perturb slighly, 
giving us "canard intervals", whose existence is stated in Theorem 1.1. 

When we consider more general case of nonconvex M and N > 1, the 
description becomes more complicated, but the main arguments still work. 
Let us assume that any two fold points lie on different vertical circles, and the 
initial point w does not lie above or below any fold point. In this case, the 
trajectory which starts in w, in forward (backward) time falls to attracting 
(repelling) segment of M, moves slowly to the right (left) until reaches the 
fold point, jumps and either leaves the base strip or falls to other attracting 
(repelling) segment of M, and the process repeats until the trajectory leaves 
the base strip (see Fig. 2). 

The main difference with the convex case here is the possibility of several 
jumps. In fact, it does not affect heuristic arguments presented above, be- 
cause they deal mostly with the after-jump rotations. However, to provide a 
rigorous proof of the main results and in particular to calculate the number 
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Figure 2: Nonconvex case: several jumps 



of limit cycles, it does not suffice to use only the ideas discussed. Instead, 
we have to perform accurate analysis of the Poincare map from T to itself, 
which is discussed in the next sections. 

3 Poincare map 

Note that the function g is bounded away from zero, so we can divide the 
system (1) by g, thus re-scaling the time: this does not change the desired 
properties of its solutions (we are interested only in phase curves), and the 
system with the new function / will satisfy the same nondegenericity assum- 
tions. Thus without loss of generality we can assume g — 1 in (1). 

Consider Poincare map P £ : T — > T. The slow motion is constant (and 
bounded away from 0), so P £ is a well-defined diffeomorphism of a circle. 
Its periodic (in particular, fixed) points correspond to closed solutions of 
the system. Denote the graph of P e by 7 e . Fixed points of the Poincare 
map correspond to the intersection points of the graph with the diagonal 
V := {y = x}. Note, that in terms of the previous section, P e (L(e)) = R(e). 

The derivative of the Poincare map in any point Xq G T can be easily cal- 
culated by integrating the equation of variations. Namely, let x = x(y; x , e) 
be the phase curve with the initial condition x(— tt;xo,e) = Xq. Then, 
P e (x ) = x(7T;x ,e) and 



dP e (x) 



dx 



(7r; x ) = exp 



1 f 



/ f x (x(y,x ,e),y 1 e)dy. (3) 



dx 



dx 
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Near the attracting (repelling) parts of the slow curve M, the function under 
the intergral sign is negative (positive), and the trajectories attract (repell) 
each other while moving in these areas. Corresponding parts of the trajec- 
tories contribute contraction (expansion) to the derivative of Poincare map. 
In "most of the cases" J f' x (x,y,e)dy 7^ and either contraction or expan- 
sion dominates, thus giving either exponentially small or exponentially big 
derivative with respect to e. 

It turns out that it is possible to replace actual trajectory in the right- 
hand side of (3) with so-called singular trajectory (or contour), which is 
defined as follows. For every point w G B \ M, which does not lie above or 
below any fold point of M, recall the description of the trajectory which pass 
through w traced in backward and forward time up to exit from B (see sec- 
tion 2). Assume that all phases of fast motion in this description are strictly 
vertical. Then we obtain a picewise-smooth curve in the base strip which 
consists of vertical segments and arcs of the slow curve M, interchanging 
each other. Call this curve singular trajectory (or contour) of w and denote 
it by Z[w). This curve is in a sense a limit (as e — > 0) of the trajectories 
with the initial condition w. The part of the contour to the right of w (which 
corresponds to the trajectory in forward time) is denoted by Z~(w), and the 
part to the left of w (which corresponds to backward time) is denoted by 



The following Lemma represents the fact that the derivative of the Poincare 
map is controlled (with given precision) by the contour of the corresponding 
trajectory. It means that the main contribution to the derivative is made by 
the segments of slow motion near the arcs of the slow curve. This contribu- 
tion dominates over the one of the jumps and the after-jump rotations. 

Lemma 3.1. Fix some 5 > 0. Fix some vertical interval J, which intersects 
attracting part of M and S -bounded from repelling part of M (and therefore 
from fold points). Let u be coordinate on J. Consider Poincare map Q : J — > 
T in forward time. Then for any w G J, 



Obviously, Z (w) does not depend on choice of w G J. The remainder term 
o(l) — > as s — > uniformly in choice of J provided 5 is fixed. 

Lemma 3.2. Let w G B\ Ms, where Ms is 5 -neighborhood of M, and y(w) 
(y-coordinate of w) is 5-far from y{G) for any fold point G. Consider actual 
trajectory that pass through w, and denote its initial condition on T by xq- 



Z+(w). 
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Then 



logP £ '(x ) 



f' x (x,y,0)dy + o(l) 



Z(w) 



(5) 



The proof of these Lemmas is given in [3] (see Lemma 5.4 and Lemma 7.3 
there), and rely heavily on the analysis of the after-jump rotations from [2] 
(see Theorem 4.3 there, which is reformulated as Theorem 4.6 in [3]). 

4 Shape of the graph of P £ 

Our goal is to describe the shape of 7 e and its dependence on e. The fol- 
lowing description goes back to Shape Lemma in [1], where it is proved for a 
particular example. 

Fix some vertical interval J + (resp. J~) in the phase space, which inter- 
sects repelling (attracting) part of the slow curve close enough to far left (far 
right) fold point and is bounded from attracting (repelling) part of the slow 
curve. (See Fig. 3.) 




— 7T 



Figure 3: Phase curves near main jump points 



Definition 4.1. A trajectory is called a duck (canard) if and only if it 
intersects J + . 
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Consider the projection of J + (resp., J~) to V along phase curves in back- 
ward (forward) time. Denote it by (resp, D~). Note that all the trajecto- 
ries that intersect are ducks. Lemma 3.1 (applied to the system with the 
time reversed if necessary) implies immediately that \Df\ = 0(exp(— Ci/e)) 
and \D~\ = 0(exp(— C^/e)) for some positive C±, G^. The trajectory with 
the inital condition x G T \ Df does not intersect J + . Therefore, it is at- 
tracted to the attracting part of the slow curve rather quickly (the speed 
is controlled by the distance between J + and the reverse jump point), and 
then moves near attracting parts of the slow curve only, accumulating con- 
traction (see the integral (3)). It also have to intersect J~ and therefore D~ . 
Lemma 3.2 implies that in this case the derivative of the Poincare map is 
exponentially small. On the other hand, appliyng the same arguments to the 
system with time reversed, we have that outside of D~, the inverse Poincare 
map P~ l has exponentially small derivative. Informally speaking, it means 
that almost all the circle in the pre-image (except for very small interval) is 
mapped into very small interval in the image, while the exceptional interval 
in the pre-image is mapped into almost all the circle in the image. 

Geometrically, this means that the graph 7 e belongs to the union II + un - 
of exponentially thin strips: vertical Il + = Df x S 1 and horizontal II - = 
S 1 x D~ . Outside of the rectangle K £ = Il + fl Il~, the slope of 7 e is either 
exponentially big or exponentially small (see Fig. 4). 





























t 





Figure 4: Graph of Poincare map 

Monotonicity arguments similar to the ones discussed in section 2 show 
that as e \ 0, rectangle K e moves from bottom-right to top- left corner, 
making infinitely many rotations. (See Monotinicity Lemmas in [1] and [2] 
for details.) 
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In this paper, we are interested only in limit cycles that correspond to 
fixed points of Poincare map (i.e. making 1 rotation along the y-axis). They 
are born (or die) when the diagonal T> tangents the graph j e . Such a tangency 
is possible only at the points where the slope of 7 £ is equal to 1. We will 
call such points neutral, applying this term both to the points on the graph 
7 e and to the corresponding values of argument (i.e. roots of the equation 
P' £ {x) = 1). Note that all neutral points belong to K £ , and therefore the 
fixed points can be born only inside K £ , thus giving us pairs of repelling and 
attracting canard cycles. (All points in K £ correspond to canard solutions 
because they lie over Df.) 

For every neutral point x, consider second derivative P"(x), and call x 
generating (resp., annihilating) if P £ (x) < (resp., P"(x) > 0). We may 
impose nondegenicity conditions, such that P £ {x) ^ in every neutral point 
x. Consider a projection A(x, y) = x — y along the diagonal V. It turns out 
(this will be discussed below) that for any given system for e small enough 
the number of neutral points is alwas the same (does not depend on e; see 
section 5) and the order of their projections under A is fixed as well (see 
section 6). In this case, actual maximal number of canard cycles is defined by 
the order of births and deaths, which is controlled by the order of generating 
and annihilating neutral points under the projection map A, and thus does 
not depends on e. This gives us the number k from Theorem 1.1. Rolle's 
Theorem implies that k < N. Neutral points depend on e continiously, 
therefore on every turn of K £ there exists an open interval of e's, on which 
the maximal number (which is 2k) of canard cycles are born. Such intervals 
accumulate to 0, and their existence is the main result of theorem 1.1. 

The rest of paper is devoted to the analysis of neutral points. We first 
prove that the number of neutral points is bounded by the number of folds 
of M and show how it can be calculated explicitly (see section 5). Then 
we discuss the order of births and deaths of canard cycles (see section 6). 
Finally, we will contstruct an open set of systems with maximal number of 
limit cycles k = N (section 7). 

5 Neutral points 

Consider a trajectory, that passes through some point w G B \ M (see the 
description in section 2). The part of the trajectory to the left from w lies 
near the repelling arcs of the slow curve; the part to the right from w lies near 
the attracting parts of the slow curve. We will say that at w the trajectory 
passes through duck (or canard) jump: the transition from unstable part 
of the slow curve to the stable one. Let S be an arc of the slow curve M 
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Figure 5: Maximal true slow curves and duck jump 



between two consequent fold points (maximal arc). It is well-known [6] that 
there exists invariant curve S £ that tends to S as e — > 0. This curve is 
called (maximal) true slow curve. It is not unique, but all such curves are 
exponentially close to each other and we can pick a suitable one. 

For every maximal arc of M, consider corresponding maximal true slow 
curve (see Fig. 5). 

Extend them in the backward time to T, and denote corresponding in- 
tersection points by u±, . . . ,u 2N (enumeration is consequent, even numbers 
correspond to the repelling curves and odd to the attracting ones; obviously, 
they should interchange). Put by definition u 2 N+i = u%- Enumerate cor- 
responding slow curves as S 1 , . . . , S 2N , and true slow curves as Si, . . . , S% N 
respectively. 

The trajectory, that passes through canard jump from S 2 i can fall after 
the jump either to S 2 i-i or to S^+i- Consider the first case. It becomes 
possible if the initial condition u belongs to the interval (u 2 i-\,u 2 i), i.e. lies 
below u 2 i. When we move u a bit upward (closer to u 2 i), the whole trajectory 
moves closer to Sf. Thus the duck jump moves to the right. It means 
that the trajectory will spend more time near the repelling part of the slow 
curve and less time near the attracting part. Therefore, it will accumulate 
more expansion and less contraction, and the derivative of the Poincare map 
increase monotonically on this interval. 2 

2 To be honest, we are cheating here a little bit: this proves monotonicity only for some 
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Similar arguments show that the derivative of the Poincare map decreases 
monotonically on interval (u2i,u 2 i+i). It follows that the Poincare map has 
picewise-monotonic derivative with exactly N intervals of growth and N 
intervals of decrease. Therefore, the equation P'{x) = 1 can have no more 
than 2N roots. This proves the estimate for the number of the neutral points, 
and therefore of the canard cycles. 

Actual number of neutral points can be calculated as the number of zeroes 
of logarithmic derivative logP e '(x). It follows from the analysis above that 
this derivative, presented as an integral (3), reaches its maximal and mini- 
mal values in the points that corresponds to the maximal true slow curves. 
Lemma 3.2 (with some modifications) implies that these values can be cal- 
culated as integrals over special contours, which contain maximal arcs of the 
slow curve. Thus the number of the neutral points equals to the number of 
sign changes of these values. It does not depend on e and can be effectively 
calculated. 

6 Order of neutral points 

Lift the rectangle K £ to the universal cover of the two-torus continuosly 
with respect to e. Pick two arbitrary neutral points £, 77 G 7 e from this 
lifted rectangle. Then we can define the difference between their projections 
A(£) — A (77), assuming that A = x — y, where x and y are coordinates on 
universal cover. (We need these precuations, because in general case the 
difference between two points on a circle is not defined.) In this section, 
we show that for any two neutral points the sign of this difference does not 
depend on e. 

The main idea is to show that the segment [£, 77] C K e is either "almost 
vertical" or "almost horizontal". In the first case, if £ is top end of the 
segment and 77 is bottom end, then A(£) — A (77) > 0. In the second case, if 
£ is left end and 77 is right end, then A(£) — A (77) > 0, and so on. 

Consider two trajectories which correspond to neutral points (call them 
neutral trajectories). Due to Lemma 3.2, they should lie near some contours 
with zero integrals (call such contours neutral as well). Note, that this implies 
that the measure of basins of limit cycles is bounded away from 0: these 
cycles lie in different areas of the phase space, which are separated by neutral 
solutions which are close to fixed neutral contours. 

Consider first forward-time parts of these contours (which are denoted 
by Z + ). They both contain the far right fold point and therefore have some 

smaller intervals of the vertical circle. Fortunately, they contains all the neutral points, so 
the proof works. See [3] for details. 
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nonempty intersection. Denote far left point of this intersection by T. To 
the left of T, the corresponding contours (and therefore actual trajectories) 
are bounded away from each other. To the right of T, the trajectories follow 
the same attracting arcs of the slow curve, and therefore attract each other. 

Consider Poincare map from some interval J' 3 T to T in forward time. 
Then the rate of the attraction is given by Lemma 3.1 and is defined by 
the integral f' x over intersection of the contours. This integral does not 
depend on e and can be calculated explicitly. The distance between these 
trajectories when they approach V is the distance between the ^/-coordinates 
of corresponding neutral points. It follows immediately that \y{C) — y(ij)\ = 
0*(exp(—C~/e)) for some C~ > 0. 

Applying the same arguments to the system with time reversed, we obtain 
similar statement for x-coordintes: \x(£)—x(rj)\ = 0*(exp(—C + /e)) for some 



Again, we may impose additional nondegenericity conditions and assume 
that C + 7^ C~ . This means that the slope is either exponentially big or 
exponentially small, what implies the necessary assertion immediately. 

This proves that the order of neutral points under the projection A is 
fixed and thefore the maximal number of canard cycles k is well-defined. It 
finishes the proof of Theorem 1.1. 

7 Duck farm 

The discussion above shows that we can translate "dynamical" questions (e.g. 
about limit cycles, Poincare map and so on) into geometrical/combinatorial 
language which involves the shape of the slow curve M and values of integrals 
of f' x over some arcs of M. As an application of this approach, we pick an 
arbitrary A > 1 and construct a system with maximal number of canard 
cycles: k = N. In fact, this example provides an open set of such systems, 
because all conditions imposed on the system during the construction of this 
example are open. To simplify the notation, we consider only the case N = 3, 
but extension of these arguments to the general case is strightforward. 

The key ingredient of the construction is the shape of the slow curve, see 
Fig. 6, top part. We demand here the depicted contours to be neutral, and 
the integrals of f x over the corresponding arcs to be equal to corresponding 
values (e.g. f HlF - f'(x, y, 0) dy = —4, and so on). 



C+ > 0. 



Consider the slope of the segment [£, r/], which is equal to 
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X 

Figure 6: The system with maximal number of ducks 



This system has 2N neutral contours, and therefore 2N neutral points on 
the graph 7 e . It follows from the previous results, that for such a system, the 
graph looks like a "staircase", where "lengths" and "heights" of the steps 
monotonically decrease (see Fig. 6, bottom part). This can be shown by 
an explicit calculation of the corresponding exponential rates that control 
"lengths" and "heights" of the steps (see the description in the previous 
section). They depend only on the integrals over the arcs which we control. 

Due to the shape of the graph of Poincare map, it follows that the order 
of the bifurcations of limit cycles is the following: first we have N births and 
then we have N deaths. During every birth a pair of cycles appear, therefore 
the number of canard cycles here is maximal and equal to 2N. Thus we have 
constructed the desired example. This proves Theorem 1.2. 
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